The role of magnetic anisotropy in tlie Kondo effect 
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In the Kondo effect, a localized magnetic moment is screened by forming a correlated 
electron system with the surrounding conduction electrons of a non-magnetic host\ Spin 
S = 1/2 Kondo systems have been investigated extensively in theory and experiments, but 
magnetic atoms often have a larger spin^. Larger spins are subject to the influence of 
magneto-crystalline anisotropy, which describes the dependence of the energy on the 
orientation of the spin due to the surrounding atomic environment^''*. Here we 
demonstrate the decisive role of magnetic anisotropy in the physics of Kondo screening. 
A scanning tunnelling microscope is used to simultaneously determine the magnitude of 
the spin, the magnetic anisotropy, and the Kondo properties of individual magnetic 
atoms on a surface. We find that a Kondo resonance emerges for large-spin atoms only 
when the magnetic anisotropy creates degenerate ground state levels that are connected 
by the spin-flip of a screening electron. The magnetic anisotropy also determines how 
the Kondo resonance evolves in a magnetic field: the resonance peak splits at rates that 



are strongly direction-dependent. These rates are well-described by the energies of the 
underlying unscreened spin states. 

A low density of magnetic impurities in a non-magnetic host metal can have dramatic effects 
on the magnetic, thermodynamic, and electrical properties of the material due to the Kondo 
effect, a many-body interaction between the metal's conduction electrons and the electron 
spin of the localized magnetic impurity^. This interaction gives rise to a narrow, pronounced 
peak in the density of states close to the Fermi energy ^ The last decade has seen a surge of 
interest in the Kondo screening of individual atomic spins, as a result of experimental 
advances in probing individual magnetic atoms by using scanning tunnelling microscopes 
(STM) ' and single-molecule transistors ' . When magnetic atoms are placed on metal 
surfaces the Kondo interaction between the localized spin and the conduction electrons is very 
strong, leading to high Kondo temperatures Tk, in the range from 40 to 200 K In order to 
probe the Kondo physics it is desirable to reduce this interaction so that the Kondo screening 
competes on equal footing with external influences such as magnetic fields. It was shown 
recently that this can be achieved by incorporating a decoupling layer between the atomic spin 
and the screening conduction electrons 

In the current study, individual Co and Ti atoms are separated from a Cu(lOO) crystal by a 

12 

monolayer of copper nitride (CU2N) as sketched in the left inset of Fig. 1 . We used a home- 
built ultra-high vacuum scanning tunnelling microscope (STM) with 0.5 K base temperature 
and magnetic fields up to 7 T to explore this system (see Methods). We probed the local 
electronic excitations by measuring the differential conductance dlldV, where / is the 
tunnelling current and Fis the sample voltage. 



A single Co atom bound to the CuiN surface exhibits a sharp zero-voltage peak in its 
conductance spectrum (Fig. 1). This peak is due to an increase in the density of states near the 
Fermi level that results from Kondo screening'. Such a Kondo peak has been observed in 
quantum dots'^'^"^ and in atoms or molecules on surfaces^'^. The peak is well described by a 
thermally broadened Fano lineshape'^ that results from interference between tunnelling into 
the Kondo resonance and tunnelling directly into the substrate'''^. For Co on CuiN the Fano 
lineshape is nearly Lorentzian, indicating that tunnelling into the Kondo resonance dominates 
over tunnelling directly into the substrate, presumably because the decoupling CuaN layer 
inhibits the tunnelling to the substrate. 

To confirm that the conductance peak is due to a Kondo resonance, we measured the spectra 
at higher temperatures (Fig. 1) and observed a rapid reduction of the Kondo peak height. 
These specfra are broadened due to the Fermi-Dirac disfribution of the tunnelling elecfrons in 
the tip, but the intrinsic temperature behaviour of the Kondo system can be determined by 
performing a deconvolution of the measured spectra'^ to remove the effect of the finite tip 
temperature (see Methods). As expected for the Kondo effect'"*, the intrinsic width of the 
resonance grows linearly with temperature at higher temperatures, but saturates at a finite 
value at low temperature (Fig. 1 right inset). The half- width at half-maximum F of the zero 
temperature peak defines the Kondo temperature Tk = F / ks, where is Boltzmann's 
constant. We find a Kondo temperature for Co on CU2N of Jk = 2.6 ± 0.2 K. 

To determine the effect of magnetic anisotropy on the Kondo screening, we applied a 
magnetic field B along each of the three high-symmetry directions of the sample (Fig. 2b and 

c). We find that the Kondo resonance splits into a double peak as in previous studies of other 
Kondo systems^'^'''''^''''''^''^. In contrast to previous studies, we observed that the peak splits 



at a rate that depends strongly on the direction of the field: when B is oriented along the 
direction in which the binding site is neighboured by two hollow sites (the 'hollow direction') 
the splitting is much less than for the perpendicular in-plane 'N direction' (Fig. 2a). The 
strong directional dependence of the Kondo peak splitting suggests the presence of large 
magneto-crystalline anisotropy^'"*. With B oriented out-of-plane the spectra are essentially 
identical to the spectra for the N direction. We note that this is true even though there is no 
symmetry between these two axes. 

We are able to determine the magnitude and the magnetic anisotropy of the local spin by 
detecting inelastic spin excitations performed by the tunnelling electrons. In addition to the 
Kondo resonance peak, the spectra show an 'outer' step in the conductance near ± 5 meV 
(Figs. 1, 2b and c). This step can be explained within the framework of Inelastic Electron 
Tunnelling Spectroscopy (lETS) as the onset of an inelastic spin excitation, which opens a 
new conductance channel ^\ The evolution of the step energy as a function of the magnitude 
and direction of an external magnetic field (Figs. 2b, c) can be used to quantify the magnetic 
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anisotropy and the spin of the Co atom . Since results obtained for B oriented in the N and 
out-of-plane directions are indistinguishable, we model the system as having uni-axial 
anisotropy along the hollow direction, which will be designated as the z axis in the following. 
We describe the single particle anisotropic spin behaviour - ignoring the many-body Kondo 
properties - with a spin Hamiltonian which is the sum of the Zeeman energy and the 
anisotropy energy"*: 

H = -gju,'B-S+DSl, (1) 
where g is the Lande g- factor and //g the Bohr magneton. The longitudinal anisotropy 
parameter Z) partly breaks the zero-field degeneracy of the eigenlevels; different magnitudes 



of m, the eigenvalue of the z-projection Sz of the spin vector S, lead to different energies 
(Fig. 3). 



We find good agreement between the measured outer steps and the corresponding excitation 
of this spin Hamiltonian (i.e. the second-lowest excitation) by using S = 3/2. The fit then 
yields g = 2.19 ± 0.09, which lies close to the free-electron g-value of 2.00, and 
D = 2.75 ± 0.05 meV, similar to the large magnetic anisotropy seen in prior studies of 
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individual magnetic atoms on surfaces ' . The uncertainties quoted are mainly due to slight 
variations in the atoms' local environments, possibly due to position-dependent strain in the 
CuaN islands'^. We note that S = 5/2 and higher half-integer values of 5* also give an adequate 
description of the measured spectra. However, all the conclusions drawn in this work about 
the Kondo and low-energy spin-excitation behaviour would remain unchanged. When an 
atom is adsorbed on a surface the spin is generally unchanged or reduced fi-om the fi-ee atom 
value. Since the spin of a fi-ee Co atom is 3/2, the spin of Co on this surface is likely 3/2. 

The sign of the anisotropy parameter D determines whether states with large or small \m\ form 
the ground state. The positive D observed here yields hard-axis (easy-plane) anisotropy, in 
which the two m = ±1/2 states have the lowest energy and are degenerate in the absence of B. 
We note that even in the case of finite transverse anisotropy the twofold zero-field degeneracy 
of these states would not be broken, due to Kramers' theorem. The observed outer step, used 
to obtain the fitting results discussed above, thus corresponds to m = ±1/2 — >■ ±3/2 transitions 
when B is applied along the z direction, see Fig. 3. 

We calculated the inelastic tunnelling spectra in the absence of Kondo screening (orange 
curves in Figs. 2b and c). The transition energies were derived by diagonalizing the spin 



Hamiltonian, equation (1), using the parameter values found above. For finite B a low-energy 
spin excitation becomes possible asm = +1/2 and m = -111 split in energy. In addition, we 
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used a scattering intensity operator that describes spin scattering by tunnelling electrons to 
model the relative step heights of the inner and outer transitions. 

We show the effect of Kondo screening on the spectra by subtracting the calculated inelastic 
tunnelling spectra fi-om the measured spectra to provide a measure of the shapes and positions 
of the split Kondo peak. This procedure results in symmetrically shaped peaks on a nearly flat 
background (Fig. 2d, e). Here we have treated the conductance as the sum of an lETS 
conductance channel and a separate Kondo channel in a manner similar to Ref 20. We stress 
that the underlying many-body physics may be more complex than the sum of these two 
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channels ' . We also note that the single-particle model of the inelastic tunnelling channel 
neglects any interference in the elastic channel due to the presence of the inelastic channel . 

The position of the split Kondo peak coincides with the calculated low-energy transition 
m = +1/2 -111 of the non-Kondo Hamiltonian of equation (1) for all measured fields and 
field orientations (Fig. 3). Here it makes no substantial difference whether we use the 
positions of the peaks in the original measurements (open circles in Fig. 3) or those that 
remain after subtracting the calculated lETS curves (crosses in Fig. 3). It is worth highlighting 
that the spin exhibits direction dependent splitting of the Kondo peak only because the spin 
{S = 3/2) is greater than 1/2. In the case of an 5* = 1/2 impurity such anisotropic behaviour due 
to crystal field effects is expected to be absent. We note that the precise nature of the splitting 
of the Kondo peak, especially for magnetic fields that are small compared to the Kondo 
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temperature, is still under theoretical debate ' ' . 



In contrast to the hard-axis anisotropy, which we find to be the case for Co on CU2N, easy- 
axis anisotropy (D < 0) would favour the m = ±3/2 doublet as the ground state of an 5* = 3/2 
impurity. In this situation Kondo screening is inhibited, as it would require Am = 3 transitions 
to be made through electron scattering. A Kondo resonance can only be formed if the 
anisotropy creates a degenerate ground state with levels connected by Am = 1 transitions (i.e. 
the flipping of a conduction electron spin). This picture agrees well with previously studied 
atomic spins Fe (S = 2) and Mn (S = 5/2) on CU2N that showed no Kondo effect. Each of these 
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was found to have easy axis anisotropy and therefore do not have a ground state doublet that 
is linked via Am = 1 spin excitations. Similar mechanisms, where crystalline anisotropy is 
responsible for creating a Kondo system from a large spin (i.e. S> 1/2), have been suggested 
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for bulk impurities ' . Recent theoretical investigations have shown that a Kondo effect can 
also occur in systems with easy-axis magnetic anisotropy, exemplified by single-molecule 
magnets, if an additional strong transverse anisotropy sufficiently mixes states with different 
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m values to create a degenerate ground state which is linked via Am = 1 transitions . 

Unlike the measurements on Co, spectra taken on individual Ti atoms (Fig. 4) show a clear 
Kondo peak but no conductance steps due to additional spin-excitations. Consequently Ti on 
CU2N can be modelled as an 5 = 1/2 system. We note that a free Ti atom has 5 = 1 in the 3d 
configuration, so in this case the binding to the surface presumably changes the atom's spin. 
Ti binds at the same location as Co on CU2N so comparable magnetic anisotropy may be 
expected. However, as shown in Fig. 4, measurements at finite magnetic fields do not show 
direction-dependent splitting. This observation confirms that a crystal field can only affect the 
Kondo resonance of an impurity that has a spin larger than 1/2. 



In summary, access to a single large-spin atomic impurity provides a new opportunity for 
studying the interplay between the Kondo effect and crystalline anisotropy. We find that a 
thorough characterization of magnetic anisotropy is essential to understanding the emergence 
of the Kondo effect. For various atoms on CuiN with spins larger than 1/2, the presence or 
absence of Kondo screening can be explained solely based on their magneto-crystalline 
anisotropy: Kondo screening can occur only if the anisotropy results in degenerate ground 
state levels connected by Am = 1 transitions. Our result will be applicable to the Kondo effect 
in other systems with large spin as well, such as magnetic atoms placed directly on a metal 
substrate or single-molecule magnets with transverse magnetic anisotropy . Further studies 
may also reveal a more quantitative understanding of the link between anisotropy and the 
strength of the Kondo screening. In addition, the ability to tune the Kondo effect by varying 
the magnitude and orientation of the magnetic anisotropy would create a new class of Kondo 
systems in which the screening could be manipulated directly through control of the local 
environment. 
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Methods 



Details of the experimental setup can be found in Ref. 1 1 . Co and Ti atoms were deposited 
onto the CuiN surface at low temperature by thermal evaporation from a metal target. These 
atoms were subsequently placed on specific binding sites by means of vertical atom 
manipulation. The differential conductance d//dFwas measured using lock-in detection with a 
50 |xV RMS modulation at 745 Hz. 

The tip was verified to have a flat density of states in the energy ranges presented, by 
observing essentially constant conductance spectra when the tip is placed over the bare CuaN 
and over the bare Cu surfaces. We determined that presence of the tip does not influence the 
Kondo system, but merely probes its density of states, by observing that the spectrum does not 
change when the junction resistance is increased. 

The temperature of the sample Tsampie was regulated using a heater, while the tip was cooled 
via strong thermal contact directly to the He refrigerator. Different refrigerator operating 
modes, each corresponding to a particular tip temperature Ttip, were used for different sample 
temperature ranges. For the deconvolution of the specfra measured when Tsampie < 1.4 K we 
used Ttip = 0.5 K, and for 1.4 K < Tsampie < 5.0 K we used Ttip = 1.8 K. With this assignment of 
tip temperatures the intrinsic peak width was found to increase without a discontinuity at 
^sample = 1.4 K. Uncertainties in Ttip up to 0.3 K were taken into account for determining the 
error bars in the right inset of Fig. 1, and hence the uncertainty in 7k. 

The precise parameter values found for the Co atoms of Fig. 2 are g = 2.16 and Z> = 2.71 meV 
for the atom with B along the hollow direction, and g = 2.22 and D = 2.79 meV for the atom 
with B along the N direction. These values were used to calculate the energy levels for the 
corresponding spins in Fig. 3 as well as the inelastic tunnelling spectra in Figs. 2b and c. 



Figure captions 

Figure 1 1 Temperature dependence of the Kondo resonance of a Co atom. The left inset 
shows a schematic drawing of a single Co atom bound on top of a Cu atom of the CuaN 
surface. Through the CuaN layer the Co atom is coupled to the electron sea of the bulk copper 
which can Kondo screen the localized spin on the atom. The main panel shows differential 
conductance (dl/dV) spectra measured with the tip positioned over a Co atom at 0.5 K and 
higher temperatures. For each spectrum the tip height was adjusted to give a 10 MQ tunnel 
junction (V=10 mV, 7=1 nA). The red curves show fits to Fano functions broadened by the 
temperature of the probing tip (Fano lineshape parameter q=n ±2). The right inset shows 
the intrinsic full width at half maximum of the peak as a function of the sample temperature T. 
The errors in these values are dominated by the uncertainties in the tip temperatures (see 
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Methods). The solid black line shows a best fit to the function [(akBT) +(2kBT]^) ] which 
describes the intrinsic width of the Kondo resonance in a Fermi liquid model' ^. For T» Jk 
the width approaches hnear behaviour with a slope a = 5.4 ± 0. 1 (red line). 

Figure 2 | Anisotropic field dependence of the Kondo resonance, a, Topographic STM 
image (10 nm x 10 nm, 10 mV, 1 nA) of four CU2N islands with single Co atoms. The two 
figures on the sides sketch the adsorption site of the two marked Co atoms (Cu and N atoms 
are depicted as yellow and green circles, respectively). Labels indicate the two non-equivalent 
in-plane directions, referring to neighbouring N atoms or hollow sites along the direction of 
the magnetic field, b, c. Black curves: differential conductance spectra taken with the tip 
positioned over the two Co atoms of panel a when a field up to 7 T was applied in the 
designated directions. Successive spectra are offset by 0.15 nA/mV for clarity. The blue 
curves in panel c show similar measurements on a Co atom where the magnetic field was 



directed perpendicular to the surface. All spectra recorded at r= 0.5 K. Orange curves: 
calculated inelastic tunnelling spectra based on the parameters obtained by fitting the 
positions of the outer conductance steps to equation 1 . d, e, Result of subtracting the orange 
curves from the black curves in panels b and c respectively, showing the change in the spectra 
due to Kondo interactions. 

Figure 3 | Energy eigenlevels for different field directions, a, Solid lines show the 
calculated energy levels based on equation 1 with magnetic field B parallel to the hollow 
direction (see Methods). Full and open circles indicate the energies of the steps and the peaks, 
respectively, of the spectra shown in Fig. 2b. The positions of the peaks in Fig. 2d are 
represented by crosses. Values are plotted relative to the calculated ground state as illustrated 
by the arrows in b. b. Same as a but with B parallel to the N direction and data taken fi-om 
Figs. 2c and e. 

Figure 4 | Kondo effect of a Ti atom: a 5* = 1/2 system. Differential conductance spectra on 
individual Ti atoms on CuiN in the absence of a magnetic field (black curve) and with an 
external field of 7 T applied along the two in-plane directions (red and orange curves, offset 
by 0.28 and 0.30 nA/mV) and oriented out-of-plane (blue curve, offset by 0.32 nA/mV). The 
7 T spectra are essentially identical because anS= 111 system cannot show magnetic 
anisotropy. All curves measured at J'= 0.5 K. 
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